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^ ' ABSTRACT. We obtained the map of the zenith brightness of the night sky in Italy, constructing a 

' simple model. The artificial sky brightness in each site is given by the integration of the contribution 

, produced by each unitary area of surface obtained by applying a propagation function to the upward 

^s^J , emission in the area. This operation is a convolution of the upward emission with the propagation 

function. In fact, the scattering from atmospheric particles and molecules of light emitted upward by 
the cities spreads the light far from the sources. In pratice, we convolved the DMSP satellite night-time 
images of the upward light emission in Italy with a propagation function, like the Treanor Law. We used 
the light emission as measured by DMSP satellite images in order to bypass errors due to differences 
in the output of cities of the same population arising when using population data to estimate upward 
flux and in order to take into account also the contribution to the sky brightness produced from cities 
Q |' outside the Italian boundaries. We chose the DMSP visible band images for their negligible number of 

I 1 saturated pixels. Italy was cloudfree as guaranteed by inspection of IR images taken at the same time. 

O ' We also evaluated the emission versus population relationship comparing the relative emissions of 

i a number of cities of various populations. The measured emissions increase quite linearly with the city 

, population in the range from 1000 to 400000 inhabitants. In a preliminary analysis, more populated 

cities seem have a lower emission per inhabitant, so that in the range from 1000 to 3000000 inhabitants 
the best fitting curve to the measured emission seems to be a power law with the power 0.8 of the city 
population. We did not find any dependence of city upward emission on the economic development of 
the area. 
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1. Introduction 

The effects of artificial lighting on the brightness of the night sky in Italy were explored 
by Bertiau, de Graeve and Treanor in 1971 (Treanor 1971, Bertiau et al. 1973), but 
since then no other systematic study followed. An inquiry done by the Commission for 
the study of light pollution of Societa Astronomica Italiana (Di Sora 1991, 1993) was 
limited to the situation of the night sky in the main Italian astronomical observatories. 
In the last quarter of century in Italy there was a great increase in the brightness of the 
night sky due to the increased number and efficiency of the lamps used for the outdoor 
artificial lighting (see e.g. Cinzano 1999b). As a consequence, awareness of the light 
pollution problem by the professional and amateur astronomical communities is greatly 



increased. A study was needed to evaluate the new situation. Mappings of sky brightness 
for extended areas were performed by Walker (1970, 1973), Albers (1998) in USA and 
Berry (1976) in Canada with some simple modelling. These authors used population 
data of cities to estimate their upward light emission and a variety of propagation laws 
for light pollution in order to compute sky brightness. Recently DMSP satellite images 
allowed direct information on the upward light emission from almost all countries around 
the World (Sullivan 1989, 1991) and were used to study the increase of this flux with 
time (Isobe 1993; Isobe & Hamamura 1999). 

In this paper we present a detailed map of artificial sky brightness in Italy. In order 
to bypass errors due to differences in the output of cities of same population arising 
when using population data to estimate upward flux and in order to take into account 
also the contribution to the sky brightness produced from cities outside the Italian 
boundaries, we constructed the map measuring directly the upward flux as detected in 
DSMP satellite night time images and we convolved it with a light pollution propagation 
law. We also studied the relation upward flux - city population in Italy and compared 
it with results obtained in other countries. 

In section 2 we describe the satellite observations. The relation upward flux - city 
population is presented and discussed in section 3. In section 4 we describe and discuss 
the method used to construct the map of artificial sky brightness. The map is pre- 
sented in section 5 together with the map of upward emission and the map of total sky 
brightness. Section 6 contains our conclusions. 

2. Observations 

We studied visual images of night time Italy obtained by the Defense Meteorological 
Satellite Program (DMSP) of the National Oceanic and Atmospheric Administration 
(NO A A). DMSP are satellites in a low altitude (830 km) sun-synchronous polar orbit 
with an orbital period of 101 minutes. Visible and infrared imagery from DMSP Oper- 
ational Linescan System (OLS) instruments monitor twice a day, one in daytime and 
one in nightime, the distribution of clouds all over the world. At night the instrument 
for visible imagery is a Photo Multiplier Tube (PMT) sensitive to radiation from 410 
nm to 990 nm (470-900 FWHM) with the highest sensitivity at 550-650 nm, where the 
most used lamps for external night-time lighting have the strongest emission: Mercury 
Vapour (545 nm and 575 nm), High Pressure Sodium (from 540 nm to 630 nm), Low 
Pressure Sodium (589 nm). The IR detector is sensitive to radiation from 10,0 fim to 
13,4 nm (10.3-12.9 FWHM). Every fraction of a second each satellite scans a ~3 km 
swath extending 3000 km in east-west direction with a resolution varying from less than 
3 km at the nadir to about 12 km at each end (Sullivan 1991). In our images each 
pixel is 2.70 ± 0.02 x 2.75 ± 0.03 km wide respectively in North-South and West-East 
directions, a mean of five by five pixels ~0.56 km wide of the high-resolution original 
images which are not distributed. The pixel values are currently relative values rather 
than absolute values because instrumental gain levels are adjusted to have a constant 
cloud reference brightness in different lighting conditions related to the solar and lunar 
illumination at the time. Due to the limited dynamic range of the satellite detectors, 
the automatic gain normally saturates the most lit pixel of the largest cities. Sensitivity 
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Fig. 1. Upward flux versus city population relationship for Italy in the range 1000 - 400000 
inhabitants. 



reaches 1CP 5 W m~ 2 sr _1 fim^ 1 (Elvidge et al. 1997). A few images were taken with 
lower gain (50-24 db) and they have only a few saturated pixels. We studied two visible 
band images chosen for the negligible number of saturated pixels and because Italy is 
almost cloudfree as guaranteed by inspection of the IR image taken at the same time. 
The images were taken the 13th January 1997 h20:27 from satellite F12 and the 11th 
March 1993 h20:37 from satellite F10 from an altitude of about 800 km. 



3. Emission vs. population relationship 

In order to evaluate the emission versus population relationship we chose a number of 
cities in cloudfree zones of the images, as seen in the corresponding IR images, and 
evaluated their emissions. To measure a city emission we summed the counts of all 
pixels pertinent to each city. So we obtained the relative emission of a number of cities 
of various populations. The population data were taken from the 1995 Italian census 
provided by Istituto Italiano di Statistica (ISTAT). 

Figure [l] show our measurements for 139 cities distributed in all Italy in the range 
from 1000 to 400000 inhabitants from our lower gain image. No sources in this plot 
has saturated pixels. The dashed line shows the best fit of u oc P. The upward emission 
increase linearly with the population in the considered range. 

In the range 1000-3000000 inhabitants shown if figures || and [| the best fit curve to 
the measured emission is a power law with the power 0.8 of the city population P. These 
results need to be confirmed as soon as many more measurements become available. 
This dependence might be due to the fact that in larger cities there are more inhabitants 
per unit area. In a district of high buildings there are less street lamps than in a same 
population district of single-family houses. This hypothesis so far was not confirmed. At 
the moment we cannot exclude that the exponent 0.8 be the result of the presence of 
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Fig. 2. Upward flux versus city population relationship for Italy in the range 1000 - 3000000 
inhabitants. 

saturated pixels near the center of more populated cities. 

A possible source of errors is that each pixel of the distributed images is the sum 
of smaller pixels and we haven't any way to check if some of them were saturated. 
This uncertainty will be solved only when original images will become available. Other 
main sources of error in our measurements are the uncertainties in determining the 
boundary of the suburbs. These errors influence mainly the measurement of the emission 
of the largest cities where the boundaries are less definite so that it is easier to include 
pixels relative to small cities in the surroundings. This causes an underestimate of the 
population relative to the measured emission. 

Figure || shows the relation between upward emission and population for Italian 
Regions when the contribution of the biggest cities is subtracted. As expected, it is 
linear. In fact, excluding bigger cities, Regions are composed of a big number of sources 
of different populations which can be considered a statistical sample. In this case the 
higher gain image was used. 

Even if the upward emission vs. city population function depends on the local lighting 
conditions, our results are in agreement with relations successifully applied or measured 
in other countries. Bertiau et al. (1973) used successfully a linear model (exponent 
equal to 1) to model light pollution in Italy, and so did Walker (1970,1973) for studies 
of California skies. Garstang adopted both models with b oc P 11 (Garstang 1986) and 
with b cx P (Garstang 1987, 1988, 1989a, 1989b, 1989c, 1991a, 1991b, 1991c, 1992, 
1993) obtaining a good fit with many observations including Walker (1977) population 
- distance data. Walker (1977) verified for a number of California cities the assumption 
of linear proportionality for street light emission only, finding a general good agreement 
with a few departures above or below the mean depending on the industrial or residential 
character of the city. Walker (1977) also measured the sky illumination produced by 
three cities of different populations, obtaining a dependence on the power of 0.8 of their 
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Fig. 3. Upward flux versus city population relationship for Italy in the range 1000 - 3000000 
inhabitants in logarithmic scale. 
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Fig. 4. Upward flux versus city population relationship for Italian Regions. 



population. Berry (1976) fitted well observations of sky brightness in city centers in 
Ontario with a propagation law for light pollution based on the approach of Treanor 
(1973) but assuming b oc P 5 . Nevertheless Garstang's linear models fit well Berry 
(1976) observations suggesting that power of 0.5 found by him be produced by extinction 
of light emitted by outskirts of large cities in propagating to the center and does not 
depend on the upward flux versus population relationship (Garstang 1989a). 

So far, we did not find any dependence of city upward emission on the development 
of the area. Cities of southern Italy have the same light output as comparable size cities 
of northern Italy, even if the former have an income per capita that is nearly half that 



of one of the latter. Bertiau et al. in the early '70 found that the city upward emission 
depended on its economic and commercial development, so they were forced to include 
in their model a development factor. There isn't any evidence in our preliminary data 
for the need of a similar coefficient now. 

4. Light pollution mapping technique 

We obtained the map of the zenith brightness of the night sky in Italy, constructing a 
simple model. If f(d) is a propagation law for light pollution giving the artificial sky 
brightness produced at a site in (x' , y') by an infinitesimal area dS = dxdy at a distance 
d with unitary upward emission per unit area, the total artificial sky brightness & at a 
site is given by: 



b(x',y') = J J e (x,y)f(V(x-x') 2 + (y-y') 2 ) dx dy (1) 

where e(x,y) is the upward emission per unit area from dS. This expression is the con- 
volution of the e(x,y) with the function f{yj (x — x') 2 + (y — y') 2 )- So our operation 
consists in convolving the satellite image giving the upward emission with the propaga- 
tion function. The scattering from atmospheric particles and molecules of light emitted 
upward by the cities spreads the light far from the sources. 

In pratice, we constructed a composite image replacing the saturated pixels in the 
higher gain image, useful to measure accurately low population sources, with the mea- 
surements coming from the lower gain image, adequately rescaled. We divided the sur- 
face of Italy in pixels with the same positions and dimensions as the satellite image. 
We assumed each area of the country defined by a pixel be a source of light pollution 
with an upward emission e X}V proportional to the counts on the composite image and 
we computed the sky brightness at the center of each pixel. In this case the expression 
(D) became: 

b ^ = XX e '^f(\J {xi ~ x h ) 2 + ( yj - y k ) 2 ) (2) 

h k 

We used for f(d) the Treanor law (Treanor 1973): 

where b is the zenith artificial sky brightness and bo is the zenith natural sky brightness 
at the site considered and produced by a city with upward emission Iq placed at a 
distance d. A, B and k are constants related to the scattering component, the direct 
beam component and the attenuation of the city light by absorption and scattering 
losses. All the constants were empirically determined by Bertiau et al. (1973) to fit in 
the best way the data of the zenith luminance due to three different cities in Italy at 
various distances from them. We checked that the ratio B /A and the coefficient k, which 
depend only on the mean conditions of the atmosphere in clean nights, were unchanged 
from the calibration of Bertiau et al. (1973). Differences between V-band and B-band 




Fig. 5. Map of the artificial sky luminance in Italy. The counts from to 100 are split into 25 
levels displayed with a multiple gray scale (bottom) in order to better show the finest details 
of the luminance distribution. Greater counts are set to 100 for clearness purposes. Contours 
of Italy are just an indication. 



propagation are under the fluctuations produced by changes in atmospheric conditions 
for clean nights (Falchi 1999). 

Using directly the upward light emission as measured in the satellite images we by- 
passed the possible errors due to differences in the output of cities of the same population 
and we take into account the contribution to the sky brightness produced from cities 
outside the Italian boundaries. This contribution was ignored by Bertiau et al. because 
bordering and near lands were sparsely populated mountain regions or underdeveloped 
countries (Albania). This approximation was far better in the early '70 than now, be- 
cause of the high increase in light emission per capita since then. A source of error should 
be the possible presence of clouds over some cities; these clouds could hide or dim the 
light received by the satellite. To avoid this we searched for an image without clouds 
over all Italian territory. We neglected curvature of the earth in our computation. This 
might produce an error for isolated areas but in strongly urbanized areas it is negligible. 
The effect of earth curvature is about 2 percent at 50 km (Garstang 1989). In computing 
the sky brightness we neglected sources outside a 200 km radius from the site. This limit 
was chosen to exclude overestimating the contribution from far cities whose effects at 
200 km are negligible due to the earth curvature. The emission coming from the same 
pixel for which we computed the artificial brightness was artificially taken away to 2.7 
km distance. This was done because we want to compute the mean brightness in the 
area, not the exact brightness in a precise site which could be sensitive to the detailed 
distribution of nearest light sources. Moreover the Treanor Law breaks down at small 
distances, causing spurious high values due to the nearest sources. To get the absolute 
values of the brightness all over Italy we need to calibrate the map with the measured 
brightness of the sky in some sites. The sites should be far from bright emission sources, 
in order to fulfil the law used. 

5. Results 

The map of artificial light pollution in Italy expressed with linear scale (artificial lumi- 
nance) is in figure ^. We also present in figure || the map of total sky brightness in Italy 
expressed with magnitude scale as Garstang (1986): 

m id = 12.603 - 2.5 log(b itj + b nat ) (4) 

where m is in V band, b is in cd/m 2 and b na t is the natural sky brightness for minimum 
solar activity assumed corresponding to 22.00 V mag/arcsec 2 . Differences with Bertiau 
ct al. (1973) due to the increase of upward emission from time and to population changes 
are under study (Falchi 1999; Cinzano & Falchi 1999). In figure [?] is shown a map of 
upward light emission as obtained directly from the satellite images. An enlargement of 
figure [s| is presented in figure |^. All the maps have been obtained from the composite 
image including low and high gain images in order to avoid saturation and to have better 
sensivity to minor sources. 

Satellite images are not calibrated so our upward flux measurements are only relative. 




Fig. 6. Map of the V-band sky brightness in Italy. The eight gray levels of the image correspond 
for 1998 to levels of <0.5, 0.5-1.0, 1.0-1.5, 1.5-2.0, 2.0-2.5, 2.5-3.0, 3.0-3.5, >3.5 mag/arcsec 2 
over the natural sky brightness. Contours of Italy are only an indication. 



7. Upward light flux in Italy. Counts go from to 171 in 16 gray levels. 
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Fig. 8. Comparisons between maps and available data. 

We tried a preliminary calibration comparing our results with measurements of sky 
brightness obtained with a specific campaign. Given the growth of light pollution with 
time (Cinzano 1999a), we used only measurements taken in the second half of 1998. 
After transforming the measurements into luminances we subtracted the natural sky lu- 
minance (see Cinzano 1999b). The comparison is presented in figure || (left). Squares are 
our measurements. Circles are available measurements respectively at Serra La Nave, 
Loiano and Ekar Observatories plotted only for comparison purposes because their al- 
titudes are higher than sea level assumed in our computations. We plan to obtain a 
better calibration after a larger sky brightness measurement campaign. 

We also tried a comparison of our upward flux map with available data of upward 
flux. Results are shown in figure || (right). Data refers respectively to Asiago, Treviso 
(from Medusa 1999, assuming public lighting 60% of total), Padova (Favero et al. 1999). 
We excluded available data from the city of Torino because it has saturated pixels in 
all our images. At the moment the resulting calibration has to be considered highly 
unconfirmed. 



6. Conclusions 

We obtained the maps of the zenith artificial sky luminance and zenith total sky bright- 
ness of the night sky in Italy convolving with a light pollution propagation function the 
upward light emission from each unitary area as measured in DMSP satellite night-time 
visible band images chosen for their negligible number of saturated pixels and for the 
cloudfree Italy. We are working to extend our mapping to other European countries. 
Figure [l^ shows a preliminary map obtained from the high gain image of 1993. Densely 
populated cities in this image have saturated pixels so the map is less accurate near 
major cities than the other maps but it could be more accurate in less populated areas. 
Better results will be presented in a forthcoming paper. 

We studied the emission versus population relationship comparing relative emission 
of a number of cities of various populations. The measured emission increases quite 



Fig. 9. Artificial sky luminance in North Italy. Enlargement of the map in Figure 5. Counts 
from to 100 are displayed in 25 levels with a linear gray scale. Greater counts are set to 100 
for clarity. 

linearly with the city population in the range from 1000 to 400000 inhabitants. More 
populated cities seems to have a lower emission per inhabitant, so that in the range 
from 1000 to 3000000 inhabitants the best fitting curve to the measured emission seems 
to be a power law with the power 0.8 of the city population. At this stage we cannot 
confirm this relation. So far we did not find any dependence of city upward emission on 
the economic development of the area. 

A comparison with the map of Bertiau et al. (1973) showing the changes in sky 
brightness from 1973 to 1997 is in preparation (Cinzano & Falchi 1999). 
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Fig. 10. Preliminary map of artificial sky luminance in parts of the European countries from 
an high gain satellite image. Due to saturation of more populated city this map is usefull only 
for a qualitative evaluation. The 16 gray levels are equispaced but the scale is not calibrated. 
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